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We derive exact results for close-packed dimers on the triangular kagome lattice (TKL), formed 
by inserting triangles into the triangles of the kagome lattice. Because the TKL is a non-bipartite 
lattice, dimer-dimer correlations are short-ranged, so that the ground state at the Rokhsar-Kivelson 
(RK) point of the corresponding quantum dimer model on the same lattice is a short-ranged spin 
liquid. Using the Pfafnan method, we derive an exact form for the free energy, and we find that 
the entropy is 3 In 2 per site, regardless of the weights of the bonds. The occupation probability of 
every bond is j in the case of equal weights on every bond. Similar to the case of lattices formed by 
corner-sharing triangles (such as the kagome and squagome lattices), we find that the dimer-dimer 
correlation function is identically zero beyond a certain (short) distance. We find in addition that 
monomers are deconfined on the TKL, indicating that there is a short-ranged spin liquid phase at 
the RK point. We also find exact results for the ground state energy of the classical Heisenberg 
model. The ground state can be ferromagnetic, ferrimagnetic, locally coplanar, or locally canted, 
depending on the couplings. From the dimer model and the classical spin model, we derive upper 
bounds on the ground state energy of the quantum Heisenberg model on the TKL. 

PACS numbers: 74.20.Mn, 75.10.Jm, 05.50.+q, 75.10.-b 



I. INTRODUCTION 

The nontrivial statistical mechanics problem of dimer 
coverings of lattices, which may be used to model, e.g., 
the adsorption of diatomic molecules onto a surface^, ex- 
perienced a renaissance with the discovery of exact map- 
pings to Ising models^. A second renaissance came with 
the search for— £ and discovery of^ a true spin liquid phase 
with deconfined spinons. In the latter case, the problem 
of classical dimer coverings of a lattice illuminates the 
physics of the corresponding quantum dimer model. At 
the Rokhsar-Kivelson (RK) point of the quantum dimer 
model, the ground states are an equal amplitude super- 
position of dimer coverings within the same topological 
sector^i^ 7 - and in fact dimer correlations at this point 
correspond to the dimer correlations of the classical dimer 
model. 

Results on classical hard core dimer models in two^ 
and higher dimensions! point to two classes of models, 
depending upon the monomer-monomer correlation func- 
tion, which is defined as the ratio of the number of con- 
figurations available with two test monomers inserted to 
the number of configurations available with no monomers 
present. On bipartite lattices (such as the square and 
honeycomb lattices), monomers are confined, with power 
law correlations^ On nonbipartite lattices (such as the 
triangular, kagome, and the triangular kagome lattice 
discussed here), monomers can be either confined or de- 
confined, and correlators exhibit exponential decay ex- 
cept at phase transtionsi 9 ' 10 ! 11 : 12 ! 13 This implies that 
while the RK point of the quantum dimer model is critical 
on bipartite lattices, so that at T=0 a (critical) spin liq- 
uid exists only at the RK point, in non-bipartite lattices, 
such as the triangular lattice and lattices made of corner- 
sharing triangles such as the kagome and squagome lat- 
tice, it has been shown that the RK point corresponds to 



FIG. 1: (Color online) A dimer covering of a portion of the 
triangular kagome lattice (TKL). The TKL can be derived 
from the triangular lattice by periodically deleting seven out 
of every sixteen lattice sites. This structure has two differ- 
ent sublattices "a" (closed circles) and "b" (open circles), 
which correspond to small trimers and large trimers, respec- 
tively. Each site has four nearest neighbors. The primitive 
unit cell contains 6 a-sites, 3 fo-sites, 6 a-a bonds, and 12 a-b 
bonds. Thick blue lines represent dimers. A typical close- 
packed dimer covering is shown. 



a disordered spin liquid. Correspondingly, it was estab- 
lished in both of these cases that there exist finite regions 
of parameter space where the ground state is a gapped 
spin liquid with deconfined spinons. Part of the interest 
in such states is the topological order that accompanies 
such ground states, and hence such states may be useful 
examples of the toric code. Interest also stems from the 
original proposals that the doped spin liquid phase leads 
to superconductivity^*^ 

In this paper, we analyze the problem of classical close- 
packed dimers on the triangular kagome lattice (TKL), 



2 



a non-bipartite lattice expected to display a spin liquid 
phase, as the first step in understanding the RK point 
of the corresponding quantum dimcr model. The TKL, 
depicted in Fig. [1] has a physical analogue in the posi- 
tions of Cu atoms in the materials CugX2(cpa)g • xH^O 
(cpa=2-carboxypentonic acid, a derivative of ascorbic 
acid; X=F,Cl,Br) 15 i 16 i 17 . We have previously studied 
Ising spins^ and XXZ/Ising spins^ on the TKL; this 
paper represents an alternative approach to the problem. 
Using the well-known Pfaffian method^ we obtain exact 
solutions of close-packed dimers on the TKL. We obtain 
an analytic form of the free energy for arbitrary bond 
weights. The entropy is | In 2 per site, independent of the 
weights of the bonds, z aa and z a b- We find the occupation 
probability of every bond is a constant j in the absence 
of an orienting potential. The system has only local cor- 
relations, in that the dimcr-dimcr correlation function 
is exactly zero beyond two lattice constants, much like 
the situation on lattices made from corner-sharing tri- 
angles such as the kagome and squagome lattices^. We 
use exact methods to find the monomer-monomer cor- 
relation function, and show that monomers are decon- 
fined on the TKL. In addition, we solve for the ground 
states of the classical Hciscnbcrg on this model. In ad- 
dition to collincar phases (ferromagnetic and ferrimag- 
netic), we find a canted ferrimagnctic phase which in- 
terpolates smoothly between the two. We obtain a vari- 
ational upper bound to the ground state energy of the 
TKL quantum Heisenberg antiferromagnet using closed- 
packed dimer picture. 

II. MODEL, THERMODYNAMIC PROPERTIES, 
AND CORRELATION FUNCTION 

In this paper we consider the close-packed dimer model 
on the TKL, a lattice which can be obtained by inserting 
triangles inside of the triangles of the kagome lattice (see 
Fig. [1]). The dimcr generating function is defined as 

z = e nv, 

dimer coverings Hj\ 

where (ij) indicates a product over nearest-neighbor 
bonds, Zij is the weight on the bond joining site i and site 
j, and riij is the number of dimers (either or 1) on bond 
ij for the dimcr covering under consideration. The term 
"close-packed" refers to the constraint that every lattice 
site must be occupied by one dimer, that is, that vacan- 
cies are not allowed. Therefore the number of sites N s i tes 
is twice the number of dimers Aimers = 12 (ij) n v- We 



allow for the possibility of different weights z a = e _/3€ ° 
for six different types of bonds a = 1,2, 3, 4, 5, 6, as de- 
picted in Fig. [21 Figure Q] shows an example of a dimcr 
covering. 

Several properties of this model, including the free 
energy, entropy, and dimer-dimer correlation function, 
can be calculated exactly using the well-known Pfaffian 




FIG. 2: Our assignment of weights z a to bonds in the TKL. 
Solid (open) circles represent a-sites (fo-sites). 

method 3 -. We begin by defining a Kastclcyn orientation 3 
(or Pfaffian orientation) for this lattice, i.e. a pattern 
of arrows laid on the bonds such that in going clockwise 
around any closed loop with an even number of bonds, 
there is an odd number of arrows pointing in the clock- 
wise direction along the bonds. For the TKL, we have 
found it necessary to double the unit cell in order to ob- 
tain a valid Kasteleyn orientation.— Such an orientation 
is shown in Fig. [3l The doubled unit cell contains 18 
sites. 

The antisymmetric weighted adjacency matrix associ- 
ated with this orientation, Aij, is a A^tes X ^Vsites square 
matrix with a "doubly Toeplitz" block structure. The 
generating function of the dimer model i s give n by the 
Pfaffian of this matrix: Z = Pf A = vdetA. In the 
infinite-size limit, this approaches an integral over the 
two-dimensional Brillouin zone: 



Ar si t cs ->oo JVsites 

where we have normalized the free energy by the tem- 
perature such that F = InZ, and where M(fe x , k v ) is the 
18x18 matrix below, 
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where, for brevity, we have written u = e tk " and v ~ e *> . 
The determinant of this matrix is independent of k x and 
k y : 

detM(k x , k y ) 

= GAzfzjz^ (ziz 4 + z 2 z 5 ) 2 {z\Za + z 3 z 6 ) 2 (z 2 z 5 + z 3 z 6 ) 2 

(4) 

Taking the logarithm and integrating over the Brillouin 
zone gives the free energy per doubled unit cell. Hence, 
the free energy per site is 
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/ = —li-i[&ziz 2 z 3 (ziz i + z 2 z 5 ) 



(5) 



x (ziz 4 + z 3 z 6 ) (z 2 z 5 + z 3 z 6 ) 




FIG. 3: (Color online) The arrows represent a Kasteleyn ori- 
entation on the triangular kagome lattice (TKL). Solid (open) 
circles represent "a" ("b") sublattices. The shaded region rep- 
resents the doubled unit cell. 



The occupation probability of each bond may be cal- 
culated by differentiating the free energy with respect 
to the weight of each bond. Let N a be the total num- 
ber of dimers on z a -bonds (as defined in Fig. [2j, av- 
eraged over all configurations of the system. Since 



z = EconfigsIL z <* a , we have N a = z a -§f-. We define 
the occupation probability of each ev-bond as p a = i g a , 
where B a is the total number of type-a bonds on the 
lattice. If TVcoiis is the number of primitive unit cells, 
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then iV si te S 

Ba = B§ = Bq = 2N ce \i s . The results, normalized by the 
number of sites in the system, are 
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Z1Z4 + Z 2 Z 5 
Z1Z4 



Z1Z4 + Z 2 Z 5 



(0) 
(7) 



Expressions for p 2 ,P3,Pb, pe follow by cyclic permutation 
of {1,2,3} simultaneously with permutation of {4,5,6}. 
The entropy can be computed by the usual Legcndrc 
transformation, S = F + /3e a N a £^ 

The behavior of the correlation functions can be de- 
duced in the same way as in Ref . S To find the dimer- 
dimer correlation functions, the standard method is to 
first calculate the "fermion" Green function, which is 
the inverse of the matrix A, Fourier-transform it to 
real space, and use the result to construct the dimer- 
dimer correlation functions. The inverse of the matrix 
A, G(k x , k y ) = [A(fc x , k y )] 1 , can be written as the ma- 
trix of cofactors of A divided by the determinant of A. 
Since det A is independent of k x and k y , the only depen- 
dence on k x and k y enters through the cofactor matrix. 
Each cofactor is at most a monomial in e lk * and e lky . 
From the rules of Fourier transformation it is easily seen 
that the real-space Green function G(x 7 y) is zero when 
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| a; | > 1 or \y\ > 1 is greater than a certain cutoff distance. 
Hence the dimer-dimer correiation function will be zero 
beyond a distance of two unit cells. This is true regardless 
of the values of the bond weights depicted in Fig. [2] This 
extremely short-ranged behavior of the correlation func- 
tion is similar to that for dimers on the kagome lattice^, 
and also to the spin-spin correlation for Ising spins in the 
frustrated parameter regime It underscores the special 
role played by kagome- like lattices (c.f. Ref . [9ll20h . 

Whereas quantum dimer models on bipartite lattices 
do not support deconfined spinons, quantum dimer mod- 
els on non-bipartite lattices can have deconfined spinons. 
The connection to classical dimer models is that at 
the RK point, correlatons in the quantum dimer model 
arc the same as the correlations of the correspond- 
ing classical dimer problem. The only non-bipartite 
lattice for which deconfined spinons have been rigor- 
ously demonstrated is the triangular lattice, by explic- 
itly calculating the classical monomer-monomer correla- 
tion function using Pfaffian methods.— On the kagome 
lattice, while no correspondingly rigorous calculation of 
the monomer-monomer correlation function has yet been 
demonstrated, there have been several indications that 
the spinons in quantum dimer models on the kagome 
lattice are deconfined (and therefore classical monomer- 
monomer correlators are similarly deconfined), from e.g., 
the energetics of static spinon configurations^, the be- 
havior of the single-hole spectral functional, and in the 
limit of easy-axis anisotropy22. We have calculated the 
monomer-monomer correlation for the kagome lattice 
dimer model using the Pfaffian approach of Fisher and 
Stephenson^ and we find that it is strictly constant, with 
M(r) = 1/4 for any r > Because the triangular 
kagome lattice dimer model maps to the kagome dimer 
model (with an extra degeneracy of 4 per unit cell), 
the monomer-monomer correlation on the TKL is also 
M(r) = 1/4 for monomers on any two &-sites, or for any 
combination of a and b sites at least three sites apart. 



III. EFFECTS OF AN ORIENTING POTENTIAL 

In the CugX2(cpa)6-a;H20 material a 15 ' 16 i 17 , the a spins 
are closer to each other than they are to the b spins, 
so the exchange couplings satisfy |J aa | > \J a b\- hi the 
classical dimer approximation described in Sec. |Vl this 
corresponds to unequal weights for dimers on ab bonds vs. 
those on aa bonds, |z aa | > \z a b\- Aside from this intrinsic 
difference in bond weights, it may also be possible to 
apply anisotropic mechanical strain to vary the lattice 
geometry (and hence the exchange couplings and dimer 
weights) in different directions. 

To obtain some insight into the behavior of the classi- 
cal dimer model under these conditions, we write z a = 
e~" e ", where (3 = 1/T is the inverse temperature and e Q 
is the potential energy for dimers on bond a. We use 
the following parametrization for the potential energy on 
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FIG. 4: (Color online) Bond occupation probabilities and en- 
tropy per site, as functions of the orienting field (3d defined in 
the text. 



each site: 



Cl — tab — 5, e 2 — £3 — e ab, 
£4 = (aa — S, 65 = 66 = € aa , 



(8) 
(9) 



where 5 is an orienting potential (i.e., an anisotropy pa- 
rameter) which favors dimers in one direction. The bond 
occupation probabilities and entropy are independent of 
the values of e a b and e QQ , and depend smoothly upon (35 
(see Fig.gjl: 



Pi = -(2 + tanh/W), 

o 

Vi = ^ (l + tanh/35) , 

P2 = P3 = 77 (4 - tanh/3<5) , 
16 

P5 = Pe = g (2 - tanh/35) , 
S 1 

s = — = — [in (64 cosh 2 (35) - 2(35 tanh (35] 

JVsitcs 18 



(10) 

(11) 
(12) 

(13) 
• 

(14) 



These results show that the TKL dimer model has nei- 
ther a deconfinement transition (as a function of e ao ~e aa ) 
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nor a Kasteleyn transition (as a function of S). It does, 
however, have a Curie-like "polarizability" with respect 
to an orienting potential. This is contrast to the situa- 
tion on the kagome lattice^ where the bond occupation 
probabilities do not depend on the orienting potential. 



IV. RESULTS FOR SYMMETRICAL CASE 



The total numbers of dimcrs on a-a bonds and on b-b 
bonds are 



_ 2 

■**a6 ^ -^dimcrs j 



(17) 
(18) 



In the absence of the orienting potential (i.e., (5 = 0), 
the expressions for the bond occupation probabilities and 
entropy become very simple: 

p a = \, a= 1,2,3,4,5,6, (15) 

s = iln2. (16) 

Note that these quantities are independent of the rel- 
ative bond weights z aa and z ao . The comparison with 
other lattices in Table [T] shows that the entropy per site 
for the TKL is the same as that for the kagome lattice. 
Although the two lattices are related, this is in fact a co- 
incidence for the following reason. The similarity can be 
seen by considering the number of 6-spins per unit cell 
which have a dimer that connects to a different unit cell. 
Because there is an odd number of sites per unit cell, this 
number must be odd, i.e. cither 1 or 3. Since the &-spins 
themselves form a kagome lattice, the same is in fact true 
of the kagome lattice. The difference is that for a given 
pattern of external dimers connecting to b spins, there 
is no further degeneracy in the kagome case, whereas for 
the TKL there are four different internal dimer patterns 
corresponding to any given pattern of external dimers 
connecting to the b spins. This means that the TKL has 
a further 4-fold degeneracy, so that the kagome entropy 
per unit cell of s ce n = In 2 becomes an entropy per unit 
cell of s ce ii = In 8 = 3 In 2 in the TKL. Since there are 9 
spins per unit cell in the TKL, this yields s — (1/3) In 2 
per site. 



where iVdimers is the total number of dimers and 
iVdimers = jJVrft*. (Of course, N aa = N± + N 5 + N 6 
and N ao — Ni + N 2 + N 3 .) Note that because there are 
twice as many a—b bonds in the lattice as there are a-a 
bonds, this implies that the dimer density is the same on 
every bond, regardless of the weights of the bonds. Since 
the number of sites is twice the number of dimers in the 
close-packed case, A^ites = 27Vdi mC rs, there are on average 
9/2 dimers per unit cell. One third of those are on the 
aa bonds, or 3/2 per unit cell. Since there are six aa 
bonds per cell, there are (3/2)/6 = 1/4 dimers per aa 
bond. A similar analysis shows that there are 1/4 dimcrs 
per ab bond. In other words, there are 1/4 dimers per 
bond, regardless of the relative weight z aa and z a b, and 
regardless of whether it is an aa or ab bond. Under the 
constraint of close-packing, the dimer densities are set 
by geometry, rather than by energetics, similar to case of 
classical dimcrs on the kagome lattic o 20 i 23 ' 24 i 25 , 

Our results for close-packed, classical dimers on the 
TKL are summarized in Table [H along with known re- 
sults for the corresponding properties on the square, hon- 
eycomb, triangular, and kagome lattices. Notice that the 
kagome and TKL are special in having simple, closed- 
form expressions for the entropies. In fact, the entropy 
per unit cell in each case is the logarithm of an inte- 
ger. On triangular lattice as well as on the two bipartite 
lattices which are shown in the table (square and honey- 
comb), the entropy is not expressible as the logarithm of 
an integer. 



Lattice 



Entropy 



Dimer corr. 



Monomer corr. 



Polarizability 



Square* 
Honeycomb^ 
TriangulapiS 
Kagome&2a 
TKL 



0.2915609 

0.161533 
0.4286 
§ln2 = 0.231049 
iln2 = 0.231049 



r 

-r/0.6014 

local 
local 



const+e 



-1/2 
-1/2 
-r/0.6014 



deconfined 
deconfined 



finite 
Kasteleyn transition 
finite 


finite 



TABLE I: Properties of close-packed dimer models on various lattices. Entropies are quoted per site. "Local" means that the 
correlation function is exactly zero beyond a certain radius - it has "finite support". The triangular, kagome, and triangular 
kagome lattices have deconfined monomers. The honeycomb dimer model not only has a finite dimer polarizability, but it has 
a Kasteleyn transition at 8 — S c . The polarizability describes the changes in bond occupation probabilities induced by an 
orienting potential 8. 
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FIG. 5: (Color online) Comparison of upper bounds on the 
ground state energy per site of the quantum Heisenberg model 
on the TKL, obtained by considering various trial wavefunc- 
tions. In the figure, we have set S = 1/2. 



The square and honeycomb lattices, being bipartite, 
admit a mapping to a solid-on-solid model^i, and there- 
fore have power-law correlations for both the dimcr-dimcr 
correlations and the monomer-monomer correlations. In 
the corresponding quantum dimcr models, these lattices 
do not support deconfined spinons. As conjectured in 
Ref. Iiol . the non-bipartite lattices have exponential (or 
faster) falloff of the dimer-dimer correlations. In the tri- 
angular, kagome, and TKL lattices, monomers are de- 
confined, which means that spinons are deconfined in the 
corresponding quantum dimcr model at the RK point. In 
fact, Mocssner and Sondhi showed that there is a finite 
region of parameter space in which a stable spin liquid 
phase is present on the triangular latticed 



V. BOUNDS ON THE GROUND STATE 
ENERGY OF THE QUANTUM HEISENBERG 
MODEL 

It is thought that the materials CugX2(cpa)6 • ajH^O 
can be described in terms of quantum S = 1/2 spins 
on the Cu atoms coupled by superexchangc interactions. 
Nearest-neighbor isotropic antifcrromagnctic couplings 
between S = 1/2 spins on a 2D lattice with sublat- 
tice structure can lead to Neel order. For example, 2- 
sublattice Neel order is favored on the square lattice, 
whereas 3-sublattice Neel order is favored on the trian- 
gular lattice.— However, on the kagome lattice and the 
TKL, quantum fluctuations are much more severe, and 
there is a possibility that they may lead to alternative 
ground states (such as valence-bond liquids). 

A valence bond state is a direct product of singlet pair 
states. Using a fcrmionic representation for the spins, 



l*W> 



(ij) 



c f c f 
V2 V lT 31 



c n c n 



|vacuum) (19) 



where riij — or 1 is the number of valence bonds on 
bond ij, just as in Eq. (fTj). 

Consider a quantum Hamiltonian with isotropic anti- 
ferromagnetic Heisenberg interactions 



H 



«7> 



Si ' Si 



(20) 



where is negative. The expectation value of this 
Hamiltonian in a valence-bond state is 

(* w | JijSi ■ Sj |* {n} ) = -^nyl^-l . (21) 

(ij) 

For close-packed dimers, the densities of valence bonds 
on aa and ab bonds are given by Eqs. (|T7|) and l|18[). 
Therefore, the total energy of the close-packed valence- 
bond "trial wavefunction" is 



Evb = -^(N aa \J aa \+N ab \J ab \) 

= — -r(\Jaa\ + 2 1 Jab |) Adimcrs 
(\Jaa\+2\J ab \)N sito s. 
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(22) 
(23) 
(24) 



This serves as an upper bound of the ground state energy 
of the quantum Heisenberg model. Of course, matrix 
elements of the Hamiltonian which connect one dimcr 
covering to another can serve to lower the actual energy 
even further. 

One may also consider a more dilute dimer state. For 
large |J aa |, one may expect dimers to preferentially oc- 
cupy a— a bonds, so that hexamers with three a — b bonds 
arc disallowed. In such a trial dimer state, the associated 
energy is 



-E'di 



lute 



(\Jaa\ + \Jab\)N s , 



(25) 



As shown in Fig. [5l this upper bound to the ground state 
energy is lower than the others for large | J aa \ . If J ao is fer- 
romagnetic and J aa is still antiferromagnetic, we expect 
another diluted dimer state, where dimers preferentially 
occupy a — a bonds, and other spins tend to be aligned 
(ferromagnetic phase). The corresponding energy is 



-Edilute+FM 



1 



\Jaa \ + g|^a&| ) N sites 



(26) 



Other bounds can be obtained by considering the clas- 
sical ground states of the Heisenberg model on the TKL 
(in which the spins are 3-vectors of magnitude S = 1/2). 
In the materials of interest, there is not yet consensus 
whether the coupling J ab is ferromagnetic or antiferro- 
magnetic. However, the Hamiltonian of the classical 
Heisenberg model is invariant under the transformation 
S b — > —S b with J ab — * —J a b, so the thermodynamics are 
independent of the sign of J a b- 

First, let us consider classical Heisenberg spins on a 
single hexamcr. By direct minimization of the energy 
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FIG. 6: (Color online) Canted state of a hexamer of classical 
Heisenberg spins on the TKL. a and (3 are the canting angles 
of the a- and 6-spins from the vertical axis. When a = /3 = 0, 
this reduces to a collinear state (which is ferromagnetic or 
antiferromagnetic depending on the sign of J a b)- When a = 
13 — 7r/2, it reduces instead to a coplanar state, in which the 
spins are all at tt/3 to each other. 




-0.2 0.0 



JaJ\Jtf>\ 

FIG. 7: (Color online) Canting angles in the ground state of 
the classical Heisenberg model on the TKL for as a function of 
coupling ratio J aa /\ J a b\- The thin line shows the canting angle 
a of the a-spins, and the thick line shows the canting angle 
/3 of the 6-spins, with respect to the collinear state, which is 
ferromagnetic or antiferromagnetic depending on the sign of 

Jab- 



of a single hexamer, we find that its classical ground 
state may be either collinear, coplanar, or canted. For 
J aa > — \J a b\/2, the ground state is collinear; the a-spins 
are aligned with each other, the 6-spins are aligned with 
each other, and the a- and 6-spins are parallel if J a b is 
ferromagnetic or antiparallcl if J a b is antiferromagnetic. 
For J aa < — | J a b |, the ground state is coplanar; the a- 
spins are at 120° to each other, the 6-spins are at 120° to 
each other, and adjacent a- and 6-spins are at 60° if J a b 
is ferromagnetic or at 120° if J a b is antiferromagnetic. 
At intermediate couplings, — | J a b| < Jaa < H^ab |/2, the 
ground state is a canted state in which neither the a- 
spins nor the 6-spins are coplanar; rather, each sublat- 



tice is canted away from Neel order, and each sublattice 
is canted away from the other. We define the canting an- 
gles of the a- and 6-spins, a and /3, such that a = [3 = 
in the collinear state (see Fig. [5]). The canting angles 
evolve continuously from 0° (collinear) to 90° (coplanar) 
as a function of the coupling ratio J aa /\Jab\ (see Fig. [TJ: 
the classical ground state has two continuous transitions. 

Now, we observe that each of these hexamer states 
can tile the kagome lattice. Therefore, the ground state 
energy of each hexamer can be used to deduce the ground 
state energy of the entire system. In the collinear regime 
{Jaa > — \Jab\/2), the collinear hexamer states lead to a 
unique global spin configuration (up to a global SU(2) 
rotation), so there is long-range ferromagnetic order (if 
Jab > 0) or ferrimagnctic order (if J a b < 0), and there 
is no macroscopic residual entropy. The ground state 
energy of the system is 

-^collinear = ^d^tml ~ 2| J a b\ )N s i tes . (27) 

6 

In the coplanar regime (J aa < — | «/a6 1 ) , there are in- 
finitely many ways to tile the TKL with coplanar hex- 
amer configurations (e.g., corresponding to 3-sublatticc 
or 9-sublattice Neel order). Furthermore, there are an 
infinite number of zero modes (rotations of a few spins 
that cost zero energy). The ground state energy is 

-E-coplanar = ~ Tj^(|>/aa| + 2| J a b\ )N s it es ■ (28) 

The physics is essentially the same as that of the classi- 
cal Heisenberg kagome model. For that model, the pre- 
vailing point of vie w 29 i 30 ' 31 i 32 is that globally coplanar 
configurations are selected at finite temperature via an 
order-by-disordcr mechanism, and the spin chiralitics de- 
velop nematic order; recently, Zhitomirsky^ 3 - has argued 
that there is an additional octupolar ordering which is, 
in fact, the true symmetry-breaking order parameter. 

The canted regime — \J a b\ < Jaa < — \Jab\/2 has the 
interesting property that in general a ^ (3, so there is a 
net magnetic moment on each hexamer. We have found 
that there are still infinitely many ways to tile the TKL, 
and that there arc still an infinite number of zero modes. 
It is possible that the zero modes cause the directions f 
the local moment to vary from place to place, destroying 
the long-range order with net magnetization; however, it 
is conceivable that the spin correlation length gradually 
increases towards infinity in going from the locally copla- 
nar state to the collinear state. The energy of the canted 
state is 

2/ 7\J aa \ 5J a 2 6 
cantod ~9lv 4 + 8|J QQ | 

- I Jab | ^(1- JlalJlb) (Jab'/JL-^ N site s ■ 

(29) 

Equations (|2T|) . (|28p. and (|2"9"|) are the exact ground 
state energies for the classical Heisenberg model on the 
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TKL. They serve as upper bounds on the ground state en- 
ergy for the quantum Heisenberg model. Figure [5] shows 
these upper bounds, plotted together with the upper 
bounds derived from dimer coverings, Eq. and (f^5j) . 
as explained earlier in this section. Notice that the up- 
per bound for the ground state energy set by considering 
dimer configurations beats the classical ground states for 
J aa large and negative (antifcrromagnctic). In this highly 
frustrated regime, we expect that the true ground state 
of the quantum Heisenberg model is significantly mod- 
ified by quantum fluctuations from that of the classical 
case. 



VI. CONCLUSIONS 

In conclusion, we have studied the close-packed dimer 
model on the triangular kagome lattice (TKL), using ex- 
act analytic methods. We find that (in the absence of 
an orienting potential) the entropy is s = ^ In 2 per site, 
regardless of the weights of the bonds, z aa and z a b- The 
occupation probability of every bond is p a = j. The 
dimcr-dimcr correlation function vanishes identically be- 
yond two lattice sites, faster than that in the triangular 
lattice, and similar to the falloff in the case of the kagome 
latticed The monomer-monomer correlation function is 
M(r) = 1/4 for r greater than two lattice constants, in- 
dicating that monomers are deconfined in this lattice. 
This implies that the Rokhsar-Kivelson point of the cor- 
responding quantum dimer model is a short-ranged, de- 



confined spin liquid. 

In addition, we find that the classical ground state of 
the Heisenberg model on the TKL is ferromagnetic (if 
J a b is ferromagnetic) or ferrimagnetic (if J ao is antifcr- 
romagnctic) when the coupling between a spins on small 
trimers is large enough compared to the coupling between 
a spins and b spins, J aa > -\J ab \/2. For J aa < -|J a6 |, 
the ground state of a single hexamer is a coplanar state, 
and the physics reduces to that of the classical Heisenberg 
kagome modcL 29 i 30 ' 31 i 32 In between, there is a canted 
classical ground state in which the a spins and b spins 
within a hexamer both cant away from the coplanar 
state. Such a state does not arise in a simple model of 
frustrated magnetism on the kagome lattice. This type 
of canted ground state of the hexamer can tile the lat- 
tice, and therefore it is the building block of the classical 
ground state of the macroscopic system. There is a cor- 
responding macroscopic degeneracy associated with the 
many ways in which this local hexamer ground state can 
tile the lattice. Each hexamer possesses a local moment; 
it is not yet clear whether the local magnetic moments 
from different hexamcrs cancel out due to the presence 
of zero modes. 
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A topological sector is defined as the following. Draw a line 
through the system, without touching any site. For a given 
topological sector, the number of dimers which cross that 
line is invariant modulo 2 under local rearrangements of 
the dimer covering. See, e.g., Ref . l34l. 

Kasteleyn's theorem may be generalized to allow complex 
phase factors in the weighted adjacency matrix: for a tran- 
sition cycle passing through sites 1,2, ...,2n, the phase 
factors must satisfy 7712)734 ■ ■ ■ »72n-i,2n = — ■ ■ ■ V2n,i- 



Complex phase factors provide a more elegant solution of 
the square lattice dimer modet^. However, they do not 
help in the case of the kagome lattice^ or TKL; we have 
found that any orientation with the periodicity of the origi- 
nal lattice violates the generalized Kasteleyn theorem, even 
if the phase factors are allowed to be arbitrary complex 
numbers. 

The expression for 5* does not simplify appreciably when 
F, etc., are substituted in. 

The long distance behavior can also be seen from the per- 
spective of the quantum dimer model of Ref. [il, in that 
each monomer removes an Ising degree of freedom, since 
it merges two hexagons. Therefore each monomer removes 
half of the configurations. 3 - 



